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■ Abstract. 



In this work we relate the curvature of distinct thermodynamic geometries to the 
response functions of any thermodynamic system with two degrees of freedom. In this 
manner it is straightforward to identify which geometry describes more accurately 
second order phase transitions. According to our results, Quevedo's metric g 11 
in general behaves better than Weinhold and Ruppeiner's, although in principle 
ambiguities might appear. It is possible to analyze the problem of describing second 
order phase transitions through the scalar curvature from a different perspective. For 
this, we propose a general criterion starting from a particular form of the curvature 
scalar. 
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1. Introduction 

Ever since Weinhold pQ [2] proposed in the seventies to use Riemannian geometry 
to study thermodynamic systems, there has been a significant number of papers 
studying the physical properties of systems by means of this formalism. This is 
done using the curvature of the Riemannian manifold representing the equilibrium 
space of thermodynamics, associating divergences of the curvature to phase transitions. 
Furthermore, Ruppeiner [3] proposed a different metric structure for the equilibrium 
space, developed from considerations about thermodynamic fluctuation theory. 
Weinhold's metric components are those of the Hessian of the internal energy function 
of the system, while Ruppeiner's are defined in terms of the Hessian of the entropy 
fundamental relation. It turns out that Ruppeiner's metric is conformally related to 
Weinhold's jl]. Recently, Quevedo [5] suggested to make use of Legendre invariant 
metrics in order to get a geometric description of thermodynamics which does not 
depend on the election of the thermodynamic potential, this approach is known as 
Geometrothermodynamics (GTD). 

Besides these thermodynamic geometries, an important and related approach in 
this context is Fisher- Rao geometry [HI [7], also called information geometry, which starts 
from statistical mechanics distributions to built the metric. A good and recent review 
on this topic is the work of Brody [7] , in which there is also an extensive bibliographical 
note and the relation between Fisher-Rao geometry and Ruppeiner's (hence Weinhold's) 
geometry is presented. We will not deal with Fisher-Rao's information geometry in this 
work, focusing our attention on the metrics arising from thermodynamics only. 

We will not describe in detail these geometries in this paper, referring the interested 
reader to the bibliography, for instance see p3, E] and references therein. We would like 
to mention only two facts about these geometric approaches: 

1. Some encouraging results have been found from testing these geometries on physical 
systems [U El El [10] , as long as some puzzling situations [HI [121 03] • In particular, 
a lot of work has been devoted to the study of black holes thermodynamic geometry 

[ia Ha us us Harm ED]. 

2. So far, every work based on the above geometric approaches has followed the same 
path, analyzing a particular system at a time. Namely, first a physical system 
and a metric are chosen and the curvature is computed, then the thermodynamic 
quantities (e.g. heat capacities) are evaluated and finally the results are compared. 

Indeed, after almost 40 years of work in this direction, still there is no agreement on 
which metric should be used for thermodynamic geometry to reproduce the physical 
properties of any system. Even more, it is not clear yet how such a metric should be 
constructed. 

In this paper we want to address the question of how to compare thermodynamic 
geometries from a general and thermodynamical point of view. In particular, what we 
will do is to obtain the curvature of all these geometries in terms of the thermodynamic 
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quantities that describe second order phase transitions. This is usually the case of black 
holes [21]. As simple as it seems, with our criterion one immediately can distinguish 
which is the best metric for reproducing second order phase transitions within the set of 
the existing ones. As a result of our work, Weinhold and Ruppeiner's metrics seem to 
be less accurate than Quevedo's metric g 11 (see [8]), as we will see in sections [2] and |3l 
However, from the form of the curvature of g 11 in terms of the response functions, some 
divergences might appear which are not related to phase transitions, although up to now 
no physical example of this situation has been found. To avoid this kind of problems, 
we suggest what should be the form of the curvature (see eq. (123]) ) and present a new 
family of metrics (eq. ( )22l) ) that describe unambiguously second order phase transitions 
as curvature singularities. 

The paper is organized as follows. In section |2] we will introduce Ruppeiner and 
Weinhold's metrics and evaluate their curvature in terms of physical quantities. In 
section [3] we will analyze Quevedo's metric g 11 and find the physical expression for its 
curvature. In section H] an ad hoc thermodynamic family of metrics for second order 
phase transitions is discussed. Finally, in section [5] we will analyze the results and 
discuss further perspectives. 



2. Ruppeiner and Weinhold's geometries. 

Weinhold's metric is defined [TJ |2] as the Hessian of the free energy U, i.e. in two 
dimensions it is the metric 

8 2 U 8 2 U 8 2 U 

First we observe that one can rewrite this metric in a more physical fashion using 
the heat capacity at constant volume, the isentropic compressibility and the isentropic 
thermal expansion (sometimes the isentropic compressibility is referred to as adiabatic 
compressibility which it is only correct for reversible processes). In fact, the heat 
capacity and the compressibility are defined as [221 [23] 

Ks = ~v[w) s = v{w^ s ' (3) 

whereas the isentropic thermal expansion is 

1/5 1 / o 2 u y 1 

as V \dTj s V \dVdSj s ' 1 j 

Hence Weinhold's metric can be rewritten as 

g w = ^rdS 2 + J—dSdV + -!—dV 2 . (5) 

C v Va s Vk s v ; 



dU 


fd 2 U" 


~dS 1 






fdV\ 


4i 


\dP) 
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Starting from this metric, one can compute the curvature as a function of the physical 
quantities Cy, Kg and as, the result being 



R 



w 



1 



V(C v ,K S ,a s ) 



(A x + A 2 a s )4 + (A 3 a 2 s + A 4 a 3 s )4 + A 5 a 3 s K S 



C 4 



v 



+ 



(A 6 a 2 s + A 7 a s )K s + A 8 a 3 s K s + A Q a s K S + A 10 a s 



c 3 



+ 



(A u a s + A 12 a s )n s + (A 13 a s + A u a s )K s + A 15 a s n s + A w a%n s 



c 2 



+ 



A 17 a 3 K% + (A 18 a s + A 19 a s )K s + A 20 a s K s C v + A 2 ia b K s 



(6) 



where 

V(C V , k s , a s ) = C v K 2 s a s [TVa 2 s - C v k s ] 2 (7) 

and Ai for i — 1, . . . , 21 are functions of the variables T, V and of the derivatives up to 
second order of Cy, k>s and as with respect to 5* and V (in particular, the A^s do not 
depend on Cy and ks). 

Usually for thermodynamic systems, second order phase transitions are signaled by 
the divergences of Cy and ks (see e.g. for ordinary systems (22] and for black holes 
[2Tj). Expression ([6]) is important because it directly relates the scalar curvature of 
Weinhold's metric to these quantities. If thermodynamic geometry wishes to reproduce 
phase transitions using the curvature, divergences of the curvature should coincide with 
divergences of the heat capacity and isentropic compressibility. Let us analyze what is 
the behaviour of the thermodynamic curvature as Cy and Ks tend to infinity. From eqs. 
(JHJ) and ([7]), we immediately notice that 

(Ax + A 2 a s )4 + (A 3 a 2 + A 4 a s )K 2 + A 5 a 3 s K S 



lim R w 

Cy— >OD 



a s n s 



T c y c 



and 



lim R 

Kg— >0O 



w 



(Ax + A 2 a s )C v + (A u a 2 s + A l2 a 3 s )C v + ^i 7 a|Cv 



a s C 4 



v 



T C ,V C 



(8) 



(9) 



where T c and V c stand for the critical values of the thermodynamic variables at the 
phase transition point. In most of the cases the limits in eqs. (EJ) and ([9]) are finite 
and show that Weinhold's geometry cannot in general reproduce the singularities in the 
heat capacity, hence under this criterion it is not a satisfactory choice, since it may 
yield to misleading results when seeking second order phase transitions. The Reissner- 
Nordstrom black hole provides an example of such a failure, since Weinhold's geometry 
is flat, contrary to the appearance of a singularity in the heat capacity [TT]. Nevertheless, 
it might happen that for some particular systems the curvature diverges at the critical 
point, but this is rather a coincidence than a general fact. 

The limit of the curvature when the isentropic thermal expansion as tends to 
infinity is of the same form as in eqs. ([8]) and (Q, although it is not clear if the 
divergence of as is related to second order phase transitions. 
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Let us now analyze Ruppeiner's metric. This metric is defined [3] as the Hessian 
of the entropy. Thus it is written in its natural variables as 

g «=-(^ + 2^ d V + §l<lV>). (10) 

In principle one cannot write this metric in terms of Cy, and as, because they are 
derivatives of U and not of S. Nevertheless, we observe that there is a relation between 
Ruppeiner and Weinhold's geometries [I], that is 

9 R =±9 W . (11) 
Eqs. (jSJ) and (ITT]) allow us to rewrite Ruppeiner's geometry as 

^-f(l dS2 + ^S dSdv + ^v dv2 )' (12) 



In Appendix Appendix A we present the curvature and its corresponding limits. There 
we argue that, as in the case of Weinhold, Ruppeiner's curvature scalar in general tends 
to finite expressions when Cy or Ks tend to infinity (the same for as). Indeed, we 
can conclude as for Weinhold's geometry, that Ruppeiner's geometry is not precise to 
describe second order phase transitions. To conclude this section, we want to point out 
that our results are obtained in general and we evaluate the limits for some diverging 
quantity keeping all the thermodynamic variables finite. This means that it might 
happen that Weinhold and Ruppeiner's metrics can actually reproduce second order 
phase transitions for particular thermodynamic systems, although this is not true in 
general [21]. In the next section we present a metric whose curvature diverges at those 
points where Cy and k$ do, making it a better candidate for a geometric description of 
second order phase transitions. 



3. Quevedo's geometries. 

Quevedo [5] has proposed a new approach in the context of thermodynamic geometry 
coined as Geometrothermodynamics (GTD). In this program metrics are derived 
imposing Legendre invariance, in order to reproduce the invariance properties of ordinary 
thermodynamics [5]. These metrics have proven to be appropriate to describe the 
thermodynamics of different systems. In particular, it is argued [8] that there is a 
metric which seems to accurately reproduce the thermodynamics of black holes or, more 
in general, of systems with second order phase transitions. Its expression in general 
coordinates, following the notation of [8], is 

where the superscript II indicates that it is a metric that reproduces second order 
transitions, $ is the thermodynamic potential, E a are the extensive variables of this 
potential, and Einstein summation convention is used. In [13] it is investigated in 
detail the problem of invariance in Geometrothermodynamics and it is pointed out 
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that, although this metric is invariant under total Legendre transformations, the same 
situation does not hold for a change of representation, so that the geometry is different 
if one chooses the energy U or the entropy S in ( fl3l) as the thermodynamic potential. 
For this reason, we now write its expression in the U representation, which reads 

$ = (S d s U + V d v U)(-d 2 s U dS 2 + d 2 v U dV 2 ), (14) 

where dxU = Using the definitions (F2J), ([2]) and the conditions of equilibrium 

dglJ = T, dyU = —P, we can rewrite this metric in terms of the physical quantities, 
obtaining 



(15) 



eJj] 



[B x + B 2 a s 1 + B 3 af) + (B 4 + Bsag 1 )*]; 1 + B 6 Kg 2 



g n ={ TS-PV){-I-^ + ^). 
From eq. f|T5l) we can compute the curvature of g]j in terms of Cy, Kg and ag, obtaining 

„ C v 

+ [(B 7 + Bsag 1 + B 9 a s 2 )K S + (B w + B n ag l ) + B^Kg 1 

+ [(S 14 + B 15 ag l )K S + Bie] Cy 1 + B 17 K S Cy 2 , (16) 

where Bi for i = 1, ... ,17 are functions of the variables T, V, P, S and of the derivatives 
up to second order of Cy, k s and ag with respect to S and V (in particular, the B^s 
do not depend on Cy and Ks). 

Assuming that the factors accompanying each of the different powers of Cy are non 
vanishing and finite at the critical points, we immediately see that 

lim R\j = ±oo, (17) 

CV— >oo 

and the same can be said about the limit when Kg tends to infinity, 

lim R\J = ±oo. (18) 

The last two equations clearly indicate that at those points where a second order phase 
transition occurs the curvature diverges. Therefore, this is the general behaviour of the 
curvature at the critical points. It is important to consider the possibility that for a 
particular system the functions Bi might be such that the limits ( TT71) and ( fl8l) are no 
longer true. However, we stress the fact that up to now every physical system that has 
been studied using g 11 is in accordance with the limits ( |T71) and (1TB]) 0, [121 El] • 

Furthermore, the curvature diverges also at those points where Cy, k s and ag 
change sign, which is also a desirable property, for these points are usually associated 
with a change in the stability of the system [2T] . 

Let us now explore the case where the potential used in eq. f|T3l) is S, i.e. the entropy 
representation. As we have already pointed out, the geometry will be different in general 
from the geometry of g\j , so we have to study it separately. In this representation, metric 
(USD reads 

gg = (U d v S + V dyS^-dljS dU 2 + d 2 v S dV 2 ). (19) 
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As in the case of Ruppeiner's metric (eq. (JT0|) ). we need an expression for this metric 
in terms of the derivatives of U in order to make Cy and k$ appear in it. However, there 
is a relation between g^ 1 and g\j (analogous to eq. (ITT]) for Ruppeiner and Weinhold's 
metrics), this relation being [13] 



" U + PV <f l (20) 
s ~ (TS - PV)T 2 

From eqs. (|15p and (T2"Uj) . g^ 1 in terms of the physical quantities reads 
„ U + PV I T Jc , 2 



9 - s ' = -— { -- dS > + — d v*y (2i) 

With this expression we evaluate the curvature and find out a result completely 
similar to eq. (JTBjl . In particular, the divergences of the heat capacity and those of the 
compressibility and their change of sign are reproduced as curvature singularities. Thus 
we conclude that g 11 is well behaved in both representations. 

For the sake of completeness, we comment that we have also examined other 
metrics which are present in the context of Geometrothermodynamics which seem to be 
appropriate to describe first order transitions [8]. For these metrics we obtained results 
for second order phase transitions which are never as good as the ones just obtained 
with g 11 , hence we will not discuss them. We conclude that, among the existing and 
already studied thermodynamic geometries, Quevedo's metric g 11 appears to be the best 
one to reproduce second order phase transitions. 

To end this section, we remark that our results are obtained in general and we 
evaluate the limits for some diverging quantity keeping all the thermodynamic variables 
finite. This means that, due to the fact that in ffT6j) the coefficients are not constant, 
there could be in principle some very unfortunate cases when Quevedo's g 11 metric 
can happen to not reproduce second order phase transitions, that would be the case 
if the coefficient of Cy in (fl6|) vanishes at the point where Cy diverges. For the same 
reason, it might occur that the curvature has some additional singularities which do 
not appear in the response functions. This means that this metric, though being the 
most accurate among the ones already existing in the literature, it is in principle not 
free of ambiguities. In the next section we will present a condition a metric should sat- 
isfy to be guaranteed that it will always reproduce the second order transitions structure. 



4. Curvature and second order phase transitions 

In this section we present a family of thermodynamic metrics which exactly reproduce, 
as singularities in the curvature, the second order phase transitions described by 
divergences of the response functions. These metrics are written in terms of Cy and Ks 
as, 

ds 2 = exp[-{aC v + bK 2 s )] {AdS 2 + BdV 2 ), (22) 
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where a and b are appropriate positive constants and A and B are arbitrary non 
vanishing constants. The curvature for these metrics reads 

R = 2eXp[ ^ +&4] {a\B{d s C v f + BC v dlC v + A(d v C v ) 2 + AC v d v C v ] 

+ b[B(d s n s ) 2 + BK S d 2 s K S + A{d V K S f + AK S d 2 v K S ] } . (23) 

Notice that the above expression has two factors which both diverge when CV and 
Ks tend to infinity (the first factor clearly diverges when Cy and ks do and the second 
factor also does, since we are dealing with vertical asymptotes of these quantities). 
Therefore, it is clear that this curvature fulfills 

lim R(Cv,k>s) — =t°° an d lim R(Cv,k s ) = ±00. (24) 

CV— >oo Kg— >oo 

Another choice to obtain (1241) is that the curvature is expressed as a linear 
combination (with positive coefficients) of the square of the heat capacity and the square 
of the isentropic compressibility, that is, 

R = Ql Cy + a 2 Kg, ai,a 2 GlR + . (25) 

One can go further and ask that the curvature reproduce also the change of signs 
in Cv and Kg- Then the equation for the curvature reads 

R = a x C v + a 2 k 2 s + a 3 + a 4 \, a x , a 2 , a 3 , a 4 G M + . (26) 

Cy Kg 

It is worth to mention again that the sole purpose of the metrics f[2"2"j) is to describe 
second order phase transitions by divergences of the curvature scalar and in principle 
not to assign a further meaning to the curvature. For instance, these metrics are not 
Legendre invariant as it is Quevedo's g 11 . It would be interesting to find out if these 
metrics have any physical origin. 

5. Conclusions and Perspectives 

In this work we calculated the general expressions for the curvature scalars of Weinhold, 
Ruppeiner and Quevedo's geometries in terms of the response functions. The limits 
when each of the response functions tends to infinity were taken for these curvatures. 
From these limits we explain in a simple and general manner why a phase transition 
described by a divergence in a response function can be identified with a curvature 
singularity. 

For Weinhold and Ruppeiner's curvatures it turns out that the scalars in general 
remain finite, which means that they cannot properly describe second order transitions. 
On the contrary, Quevedo's curvature diverges at those limits, thus indicating that it is 
a better candidate to describe second order phase transitions. 

Moreover, we found a family of metrics fl22|) whose curvature (123]) readily fulfills 
the task of diverging if and only if the response functions diverge. With the use of 
these metrics, it is possible to describe unambiguously second order phase transitions 
as curvature singularities. 
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We would like to stress that we are only focusing on the behaviour of 
thermodynamic geometries with respect to second order phase transitions points, leaving 
aside any other physical or invariance property of the metric. 

We believe that this work could be the starting point to a more detailed and wider 
analysis of phase transitions from the geometric point of view, that eventually will lead 
us to a geometric definition of phase transitions. 



Acknowledgments 

The authors want to thank Prof. H. Quevedo for insightful discussions. A. B. would like 
to thank ICRANet for financial support. F. N. acknowledges the support from DGAPA- 
UNAM (postdoctoral fellowship). This work was supported by DGAPA-UNAM grant 
No. IN106110 and CONACYT grant No. 166391. 



Appendix A. Curvature expressions and limits 

The curvature of Ruppeiner's metric ( 1121) is 



R R 



•1 



V(C v ,K S ,a s ) 



(Ax + A 2 a s )4 + (A 3 a 2 s + A A a s )K 2 s + A 5 a s K S 



Ct 



v 



+ 



(A 6 a 2 s + A 7 <x 6 s )k% + A%a s K 2 s + A 9 a s K S + A w a s 



cl 



+ 



{Anal + A l2 a%)4 + {A n a% + A u a s )4 



+ (Ai 5 <4 + ^16<4) K I + A n a s K s 



r 2 



+ 



Aiza%K% + (A 19 a s + A 20 c4)4 + A 2 ia b s n 2 s 



where 



V(C V , k s , a s ) = C v K 2 s a s [TVa 2 s - C v k s \ 



(A.l) 



(A.2) 



and Ai for i = 1, ... ,22 are functions of the variables T, V and of the derivatives up to 
second order of Cy, k-s an d as with respect to S and V (in particular, the AiS do not 
depend on Cy and ks)- Notice that here the functions Ai are different from the ones in 
OH]). The limits when Cy and Ks tend to infinity are of the same form as for eqs. (1S1) 
and (Q. Therefore, Ruppeiner's approach suffers the same inaccuracy as Weinhold's. 
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